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We study effective-medium approximations for linear composite media by means of a path integral 
formalism with replicas. We show how to recover the Bruggeman and Hori-Yonezawa effective- 
medium formulas. Using a replica-coupling ansatz, these formulas are extended into new ones which 
have the same percolation thresholds as that of the Bethe lattice and Potts model of percolation, 
and critical exponents s = and t = 2 in any space dimension d > 2. Like the Bruggeman and 
Hori-Yonezawa formulas, the new formulas are exact to second order in the weak-contrast and dilute 
limits. The dimensional range of validity of the four effective-medium formulas is discussed, and it 
is argued that the new ones are of better relevance than the classical ones in dimensions d = 3, 4 for 
systems obeying the Nodes-Links-Blobs picture, such as random-resistor networks. 



PACS numbers: 05.10.-a, 05.40.-a, 05.50.-)-q, 72.70.+m 



I. INTRODUCTION 



Among various effective-medium formulas used to model the effective behavior of random conducting linear com- 
posites, the symmetrical Bruggeman formula is undoubtly the most popular. Applied to an insulator/conductor 
binary mixture, it predicts a percolation-like transition |^-^ for a volumic fraction of conductor pc = 1/d, where d 
is the space dimension. The critical exponents are s = t = 1, and its critical properties are thoroughly discussed in 
Ref. This formula can be interpreted in two ways. On one hand, Milton has shown that it yields the exact 
effective conductivity of an ad hoc ideal medium built with a particular hierarchical structure. On the other hand, the 
Bruggeman formula can be seen as a first (one-body) self-consistent approximation to general disordered symmetric 
cell-materials |^, to which systematic corrections could be worked out. However, the Bruggeman approximation is 
very different from a mean-field theory of random conducting media. Indeed, an exact mean-field calculation on the 
Bethe lattice predicts a percolation threshold Pc ~ l/(2d) and exponents 5 = and t — 3. These exponents 

are exact for c? > 6, as well as the asymptotic behavior of the threshold when c? — > oo (at least, for the hypercubic 
lattice 1^] to which a continuum theory naturally compares 0). These values are also obtained in a more systematic 
mean- field theory for random resistor networks [ p^ . The remarkable discrepancy between the mean- field results and 
Bruggeman's formula indicates the ambiguous status of the Bruggeman theory. As a matter of fact, in spite of various 
(mostly perturbative) investigations |l^-p^ in order to precise its theoretical status, the reasons for the peculiar 
critical behavior of Bruggeman's formula are not completely cleared up. More surprisingly, another self-consistent 
effective-medium approximation [ p^p6| due to Hori and Yonezawa (HY), obtained for the same type of media by 
means of a completely different approximation scheme (and later derived by functional methods ^^), exhibits the 
same exponents s — t = 1 and a similar threshold behavior pc = I — cxp(— ~ 1/d. 

Apart from phenomelogical variants, and up to our knowledge, the Bruggeman and HY effective-medium formulas 
are the only ones obtained from the equations of electrostatics in continuous media which are able to describe, at least 
qualitatively, the overall features of a percolation transition in any dimension. One intriguing question concerns the 
possibility of deriving alternative effective-medium formulas from a continuum formulation, which do not lead to the 
seemingly unavoidable values s = t = 1 and Pc ^ 1 j d. As we show in this paper, such a possibility exists. Our starting 
point is the path integral approach recently put forward by Barthelemy and Orland , where the effective-medium 
problem is recast in a functional form. The problem reduces to compute a free-energy: roughly, the logarithm, 
averaged over the disorder, of a functional integral of Boltzmann-like weights, over allowed field configurations (which 
include boundary conditions). The average of the logarithm is carried out with the replica method (already used 
in Rcf. jl^). In Ref. |l^, the authors showed that the path integral formulation allows one to easily recover the 
second-order weak-disorder expansion of the effective permittivity of nonlinear composites |l9] . 
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However, this formulation has not yet been used to derive self-consistent estimates. In this paper, we show how 
this can be done. After a presentation of the functional approach to the homogeneizati on p roblem, and of the 
replica method (Sec. ||), we discuss self-consistent effective-medium approximations (Sec. |ni| ). As usual in such 
approximations, a background reference medium is introduced under t he for m of an ansatz for the energy of the 
system, whose parameters are 



SeTll 



The new feature here is that the 



to be determined self-consistently (Sec. 
ansatz contains a replica-coupling term, whose significance is explained (Sec. |1II B ). The self-consistency con dition s 
to determine its parameters are next discussed, and two types of effective- medium formulas are identified (Sec. IIIC): 
one in which the replica couplings are cancelled (hereafter referred to as "type 1"), and the other one with non-zero 
replica couphngs ("type 2"). Two different approximations are then worked out for each type (Sec. IV). It is found 
that type 1 generates the Bruggeman and HY formulas, whereas type 2 brings in two new effective-medium formulas 
which are "replica coupling counterparts" of the previous ones. They possess exponents s = 0, t = 2, and a threshold 
Pc ~ l/(2rf) (Sec. I V B 3| ) . These new formulas are discussed in Sec. 0, where numerical results are presented before 
we conclude in Sec. VI. 



II. PATH INTEGRAL FORMULATION OF THE PROBLEM 



The effective properties of a random conducting medium can be defined with the help of the total dissipated power 
in the medium po|-p2|. In terms of the electric field E{x), the dissipated power w in the system of volume V reads 

W[E]= [ dxw^{E{x)), (1) 
Jv 

where Wx is the local power density. 

Hereafter, we take the volume V of the sample equal to one. In heterogeneous materials, Wx depends on constitutive 
parameters randomly varying from point to point. For linear conducting media with j{x) = (t{x)E{x), where a is the 
local random conductivity and j is the electric current, we have 

w,{E{x)) ^aix)E^ix)/2. (2) 

In the analogous effective permittivity problem the dissipated power is replaced by the stored energy e{x)E'^{x)/2 (e 
is the permittivity). For this reason, we shall abusively refer to Wx as the "energy density" hereafter. In the nonlinear 
problem, Wx{E) is a non-quadratic function of E. 

An alternative to solving Maxwell's equations is to minimize the total energy W subjected to the two constraints 
|^,|l|: (i) E = —Vcj) and (ii) E — Eq; here, the bar stands for a spatial average, and Eq is a constant applied electric 
field. The minimum, W*[Eo^, is expected to be self-averaging, as occurs for the free-energy in disordered systems. 
We can therefore write 

W*{Ea) = ( min W[E]) (3) 



E = Eo 
i B=-V0 



where the brackets (•) denote the disorder average. W*{Eo) is the energy in a homogeneous medium characterized 
by an effective constitutive law Q 

dW*{Eo) 

{]) = — — = CTcff-Eo. (4) 

The second equality defines the effective conductivity of the medium. 

The problem thus reduces to computing the average of the constrained minimum of a functional of the electric field. 
The electric field derives from a potential and has a fixed mean value. We can rewrite the constrained minimum in 
(^ using a path integral 



min W[E] = - lim ^In [ VEV(j>S{E + \7(p)6{E - Eo)e-^^^^^l 

e=Eq 0— too p J 



E^ 
E = - 



(5) 



The minimum can be interpreted as the ground state energy associated to the partition function 

Z = [ VEe-^'^^''\ (6) 
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where we have used the shorthand notation 



VE = VE S{E - Eo) / V(P 6{E + V0) (7) 



for the constrained functional measure. We need to compute the average of the logarithm of (^. In order to proceed, 
we introduce rephcas p3| , p4[ and use the identity (InZ) = lim„^o((-^") — 1)/", hence 

= - hm hm - 1). (8) 

The limits do not commute. The equivalent form 

W* = - lim lim 4t In (Z") (9) 

can be used as well. The replica method relies on the fact that one can easily compute the replicated partition function 
(Z") for n integer, and subsequently take the limit n ^ 0. The main quantity of interest therefore is 

/n 
J]P£;" ^e-^EL,^[^°l). (10) 

Q = l 

Denoting the replicated measure by 'D{E"^^ = IIq^i T^E", the average (Z") can be written in terms of an "effective 
Hamiltonian" 

(Z") = y e~'^^% (11) 

with 

W. = -iln(e"''i::^.^l^°l). (12) 

For simplicity, we restrict ourselves to cell materials where the local properties are statistically uncorrelated from 
site to site. Volume integrals may then to be identified with sums over sites (each pertaining to one cell) according to 
the correspondence J dx v , where v is an infinitesimal cell volume (which defines the microscopic correlation 
length of the problem). Then, TCe simplifies to 

H, = -i|^ln^e-^"E.-(^°(^))y (13) 

Note that our discussion in Sec. Ill will be specialized to binary disorder for which the constitutive parameters can 
take only two values (but the proofs are general). That is, we assume that the local energy density is distributed 
according to the probability distribution 

P(w^w^{E)) =pd{w-wi{E)) +qd{w -W2{E)). (14) 

(where q = 1 ~ p). With this choice. 



(15) 



The above formalism applies to any form of the energy density, and in particular to nonlinear media p5|-p8[. A 
method for extracting from the path integral the second-order weak-contrast perturbation expansion of the effective 
potential W*{Eo), for nonlinear media, has been introduced in Rcf. p8|. 



III. PRINCIPLE OF SELF-CONSISTENT APPROXIMATIONS 



In this paper, we consider the linear problem only. This section is devoted to self-consistent approximations to 
W*. We first present the principle for building such approximations through the introduction of a trial Hamiltonian. 
Then, we discuss the choice of a trial Hamiltonian with replica couplings. Finally, we explain how to exploit these 
replica couplings in order to obtain two kinds of self-consistent formulas. 
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A. Overview 



The common ingredient to the approximations discussed below is the introduction of a hncar comparison medium 
described by a trial Hamiltonian Tig which is quadratic in the electric field and non-random, e.g. the one-parameter 
ansatz 

no = ^ fdxY^E'^^x), (16) 

•' a 

where do > is to be determined by an appropriate self-consistency condition. This Hamiltonian is that of a 
(replicated) homogeneous medium, but without couplings between replicas. Its meaning and that of other possible 
choices with replica couplings are discussed below. 
The partition function (Z") can be rewritten as 

or, with another notation 

(Z") = (e-^(«'-«''))^Zo, (18) 

where Zq is the partition function associated to Tio, and (•)q stands for the functional average with weights e~/^'^° /Zq. 
Equ. ^ thus reads 

W* = Wo + AW, (19) 

where 

WoiEo) = -lun^lnZo, (20) 

AVK(£;o) = - lim — ln/e-'^(^=-^°)\ . (21) 
nB \ / 

The quantity /SW{Eq) is difficult to compute (an exact evaluation would lead to the exact result for the effective 
conductivity), and we have to resort to approximations. 
A natural self-consistency condition for TLq is 

AW{Eo) = 0, (22) 



which completely determines Tip in the case where it depends on one single parameter, as in (^^. For more general 
choices of Tip with several free parameters, ( p2|) only provides a relation between these parameters, and additional 
considerations are in order to determine them all. First of all, we have to precise the form of the ansatz to be used 
in our calculations. 



B. Replica couplings and choice of the ansatz 



We deduce here the form of the trial Hamiltonian TYq from an analysis of the effective Hamiltonian Tig. Eq. (|lz 
shows that the effective Hamiltonian is non-random but that the average over disorder introduced a coupling between 
different replicas. The meaning of these couplings is more transparent if we carry out an expansion of ( |l2| ) around 
the average field E = Eq as in the weak-contrast expansion lia]. With di = d/dEi and Ai?" = i?" — Eq we have 



He = n{wx{EQ)) 
1 '" 



E 



dx a,^AE°'{x)AE°'{x) 



dxdy c-f {x 



-y)AE^{x)AEUy) 



(23) 



where 
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(x-y)^ {d^Wx{Eo)djWy{Eo)) - {diWx{Eo)){djWy{Eo)). 



(24) 
(25) 



The first non-zero replica-coupling term is proportional to /3c^^^ . We thus see that the coupling between replicas acts 
only within clusters defined by n-point connected correlation functions c^"^ , and accounts for the fluctuations of the 
electric field in these clusters. The replica coupling would vanish if there were no disorder at all. In the limit where 
the size of the region defined by c*^^) shrinks to zero - which means that the system is observed at a macroscopic level, 
we can approximate 



(2) 

y) ~ vc\^'{0)5{x - y), 



and we recover the expansion 



He 



n{wx{Eo)) + i 



dx 



Ai?r(x)Ai?"(x) - «/3 Vc^f (0)Ai?r(x)Ai?7(a;) 



a, 7 



(26) 



(27) 



which could directly be obtained from (|13|). The presence of v in front of the replica coupling term is the macroscopic 
remnant of a microscopic average having been taken within a two-particle cluster, of center x and volume v. This 
discussion therefore enlightens a relation between replica coupling and the electric field fluctuations within clusters. 
Expansion (E^) suggests a two-parameter replica-symmetric ansatz of the form 



where 



Ho — 



2 ^ 

a7 



dxM'^^ETEl 



(28) 



(29) 



The free parameters are ctq and Q. Note that Q has the dimension of a squared conductivity, because it is related to 
a quantity relative to two points. For simplicity, the ansatz M is diagonal in the euclidean vector space. However, 

we tried calculations with a tensorial structure reproducing that of a^j and Cj-^"* in Eqs. (24), ( |25|) ; but, apart from a 
different normalization for Q, no differences showed up in the final effective- medium theories (as far as linear media 
are concerned). 

An interesting feature of the ansatz ( p9| ) is that, though being non-random, it embodies underlying disorder through 
its replica couplings. In order to understand this point, we compute Wq{Eq) given by (|20| ) 



W^{Eo) 



71/3 J 



After writing E = Eq — V(/), and going to the Fourier transform of i 



WoiEo) 



— — lim — In 



(Det M) 



^, we arrive at 



(30) 



1 



a^{l-Q/al)E^. 



Carrying out the derivative of (pO[) with respect to ctq, we obtain 



El 



^0 



(31) 



(32) 



where volume averages E'^ have been replaced by statistical ones, the microscopic size f ^^"^ being much smaller than 
that of the system, y^/'^ = 1. All the replicas are equivalent, and the functional average (■)o selects in the limit 
(3 ^ oo the real field in the medium. Hence, setting AE — E — Eq, the previous equation leads to 



(AE^) Q 



El 



(33) 



which implies that Q > 0. When Q ^ 0, the electric field fiuctuates in the medium, whereas it is uniform when 
(5 = 0. The ansatz Tio therefore represents a medium which is homogeneized (because it is non-random), but which 
nonetheless accounts for field fiuctuations. We thus expect new effective medium approximations when the replica 
coupling Q is non zero. 
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C. Self-consistency 



Up to this point the discussion focused on the ansatz itself, without referring to Tig. In particular, ctq was treated 
as a mere number. We now discuss what happens when self-consistency is imposed, within some approximation 
scheme. The medium is made of N phases labelled by v, of respective conductivities ct^ and volume concentrations 
p^. The self-consistency relation AW{Eo) — which imposes constraints on the ansatz, determines Q as a function 
Q = Q{ao, {a^}). Then W* = Wq and, with 



Coff = CTQ 



(34) 



Suppose now that ctq — aQ{{a^}) is determined by an additional condition (to be precised below). Using the exact 
formula |3^] (cf. Appendix |^ 



^0 



(35) 



the fluctuations of the electric field deduced from (p4) can be written 



{AE^ 



f^off(c^o) 




1 +^-- 



dQ ^ f dQ 

9(70 ^ ^ ' 



(36) 



where the last derivative is performed at constant CTq. This expression distinguishes between different contributions to 
the field fluctuations: (i) the first term represents fluctuations coming from the "macroscopic" background effective 
medium uo; (ii) the second one is that already found in (^), and would be the only one if Q were independent from 
(TO) and if (Tq were equal to (cr), the trivial value corresponding to a non- fluctuating reference medium for a multiphase 
composite, cf. Appendix (iii) and finally a third term comes from the dependence of Q on ctq and a^. Both last 
terms are, according to the interpretation of replica coupling developped in the previous section, of "microscopic" 
origin. 

We now turn to the determination of cro({cri/}). A first obvious self-consistency condition for is Q = 0, so that 
(Teff = ctq. The effective-medium formulas obtained this way are referred to as "type 1" hereafter. As is shown below, 
to this type pertain the Bruggeman and HY formulas. 

"Type 2" effective-medium formulas are obtained by taking cto as the solution of a'^^{ao) — 0, and by using this 
value in GcS- According to (|3^), this procedure makes the effective-medium insensitive to the fluctuations generated 
in the reference medium ctoi so that relevant fluctuations only come from Q. 



IV. TWO APPROXIMATIONS 



In this section, the ideas introduced above are used within two different approximations to AW{Eq), based on the 
ansatz (|8|), (^). For each approximation to AW , "type 1" and "type 2" formulas are obtained. Herafter, q = Qju^^ 
so that 

o-off = o-o(l - g). (37) 



A. One-impurity approximation 

We first consider a "one-impurity" (or "local") calculation. The Bruggeman formula emerges as the "type 1" 
effective-medium formula in this approximation, which is not suprising since it can be seen as a one-site (self-consistent) 
theory |l|,|l|. 
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1. Approximation scheme 



The approximation for AVl^(£^o) [Eq. (pT])], detailed in Appendix is a one-impurity approximation where inter- 
actions between different points are ignored. Let us denote by wq the trial Hamiltonian density, which depends on all 
the rephcas, defined from (pSf), (|29|) by 



Ho = y dxwQ[E{x)]. (38) 
Here and in Appendix the notation [•] indicates a dependence with respect to all the replicas. Setting 

Au;,[^(a;)] =^ (a:)) ~w;o[^(a;)], (39) 

a 

the onc-impurity approximation results in 

Because of statistical translation invariance, the final result is independant of the point x. The right-hand side can 
be computed exactly for any potential in the limit [3 ^ oo using a saddle-point method. Setting Act = (t — ctq and 

M-fl + ^V. (41) 



dao, 
we arrive at 

^W{E^) = ]^{l^a^i)El + 'i{a^l)El (42) 

The condition /^W{Eo) = yields 



2. Type 1 formula: Bruggeman's 



Letting q = amounts to imposing (Aafi) — 0, which is nothing but the Bruggeman equation 

CT - (To 



a + {d~ l)(To 



= 0. (44) 



The Bruggeman equation can also be written (fi) = 1, or cto = (f/i) if d 7^ f , or cto = (afi)/ {11). The last expression is 
suitable for computing ctq iteratively (starting, e.g., from (To — (o")) in ^^ly dimension. The Bruggeman conductivity 
(Teflf = (To possesses a percolation threshold Pc — ^/d, and critical exponents s ^ t = 1 [Q. 
The fluctuations computed from ( |35| ) read 

3. Type 2 formula 

We now let g ^ and given by ( p3[ ) and (Toff((To) = (To(1 — q{ao))- The equation a'gg{ao) — reads 

(To^M^ + Mm^J^i!)MV (46) 
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Like Bruggeman's, this equation is easily solved by iterations starting from ctq = (c). The iterations then always 
converge to the physical solution, which we denote by ctq. The effective conductivity thus is CToft = o'o(l ~ ^('''o))- '^'^ 
study its critical behavior, we consider a binary mixture, where a — ai with probability (1 — p), and a = a2 with 
probability p. In the conductor/superconductor limit where f72 — > oo we find, setting pc = l/{2d — 1), 



pid-1) 



1 - P/Pc 



- 1 



{P < Pc) 



(47) 



and 



I -dp- V(l-P)(l-P/Pc) 



1)2 



(P < Pc)- 



(48) 



The critical concentration pc can be interpreted as a percolation threshold, and is the same as that obtained in the 
mean- field model on a Bethe lattice with connectivity z — 2d. Since Ueff = 2{2d — l)ai/{d — 1) ~ {Pc — p)'^ for 
P ^ Pc, the superconductivity exponent is s = 0. Note however that ctq displays a square-root cusp at p = Pc- The 
critical behavior for p > is obtained by exmining the insulator/conductor mixture where a2 is finite and CTi = 0. 
Then 



Ctq = 0-2 



CTcff = ^2 



P/Pc - 1 
2(d-l)' 

{p/Pc - I)' 
4(d- l)2p 



{P > Pc)- 



(49) 
(50) 



Since tToff (p — Pc)^ for p Pc, the conductivity exponent is t = 2. 

For the special case oi d = 1, fi — <To/a so that (|4^) reduces to do = a-^d q 0. Therefore, dcfr = 

which is the exact result. Like Bruggeman's, the new formula is also exact to second order in the contrast, in any 
dimension 



o-off = (cr) 

and in the dilute limit where (e.g.) P2 ^ 1 



1 - 



fcff = CTl 



1+d- 



d{a)^ 



<J2 — O"! 



a2 + {d - 1)cti 



(51) 



(52) 



In the discussion (Section it is argued that because of its exponents s 7^ t, and because it is less trivial than 
the Bruggeman formula (especially in the insulator/conductor case where it does not reduce to a straight line), this 
formula may constitute an easy-to-handle alternative to the latter in dimensions rf > 3. Graphical comparisons 
between different effective-medium formulas are discussed in Sec. 0. 

B. Cumulant series approximation 

In this section, we show how to recover by means of a cumulant approximation to /SW the effective-medium formula 
of HY, together with its "type 2" counterpart. 



1. Approximation scheme 
We consider the first-order cumulant approximation 

We have then 



(53) 
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AW{E^) ~ lim - {He - no)o 



(54) 



As is shown in Appendix C, the calculations here involve an expansion in a series of the cumulants of the distribution 
of (T, whose significance has been discussed at length in the original paper by HY After some algebra, we obtain 
(cf. Appendix C) 



AW{Eo) = [1 + dho{l/{dao))] + dqho{l/{dao))}Ei 



where 



ho{z)= / dMe-"ln(e-"'^^). 

The family of functions hk is defined in Appendix C. The self-consistency AW{Eo) = now yields 

1 

1 



dho{l/idao)) 



(55) 



(56) 



(57) 



2. Type 1 formula: the Hori- Yonezawa formula 
We first consider the case with no couplings between replicas, i.e. q ~ 0. The HY formula for ctq reads 

/io(l/(dao)) =-i 



(58) 



and the effective conductivity is (JcS ^ (Jo- It can be show n tha t (Toff displays a percolation threshold pc = 1— exp(— l/d), 
and exponents s — t — 1 . Applying (|35| ) and using ( C12|) , the fluctuations read 



(AE^) _ 2 + dhi{l/{dao)) 
El " 'l + dhi {l/{dao)y 



(59) 



3. Type 2 formula 

We now consider q ^ and determined as a function of o-q by (^^. Then 

1 



O'cff(o'o) = fO 



dho{l/idao)) 



The equation for do is o'^ff(o'o) ~ 0; that is, with ( C12D 



l /ii(l/(dao)) ^ 

d/l2(l/(dao)) ■ 



(60) 



(61) 



In order to study the critical behavior of (TcS, we consider again a binary mixture where a = ai with probability 
(1— j>), and (T = (72 with probability p. In the conductor/superconductor case where (72 — > cx), we have /io(l/((i(7o))) = 
ln(l — p) — (7i/((iCTo) and a similar equation for hi, so that ( |6l| ) reduces to a second-degree polynomial equation. Its 
physical solution reads 



2ai 



y/l + 2d\n{l -p) 1 + ^l + 2dln(l-p) 
It is defined for p less than a critical value 



(62) 
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-l/(2d) 



(63) 



This percolation threshold is the same as the one obtained in the Potts model at the mean-field level, and in the 
mean-field theory of Ref. |12|. Reporting (^) into (60), we arrive at 



4(7l 



1 + v/l + 2dln(l-p) 



(P < Pc)- 



(64) 



Since a^g oc (pc — p)" for p Pc, the superconductivity exponent is s = 0. 

In the opposite insulator/conductor case, where (Ti = and 1J2 is finite, the solution for p > pc can only be found 
perturbatively around the percolation threshold. Expanding the logarithm in ho and hi as 



In 



(1 -p)+pe-"'^^/('^'"°) ^\n{l~p)+J2 

l>l 



(-1) 



1-p 



^-lurT2 / (da-o) 



and defining 



we find that 



and that 



CT2 



-{p - Pc) + O {{p - pcf) 



CToft 



" A<PA{pJ{l-p,)) 

-{P-Pcf + 0{{p-pcf) {p>Pc), 



CT2 



4d2v4(p,/(l-p,)) 



(65) 



(66) 



(67) 



(68) 



where the conductivity exponent is t = 2. Hence, as in the previous "one-impurity" approximation, the replica- 
coupling ansatz yields critical exponents s = 0, t = 2, and an asymptotic dependence of the threshold Pc ~ 1/(2^) 
when 2d ^ 1. One can easily check that this new "type 2" effective- medium formula is exact to second order in the 
weak-contrast limit and in the dilute limit. 

After a few manipulations, we now obtain with (|35|) and (|57|) 



(69) 



V. DISCUSSION 



We plot in Fig. (resp. ^ the "type 2" scaled conductivities (Tcs/<^i versus p, the volume fraction of material 2, for 
a dielectric ratio a-ijai — 10 (resp. 02! o\ — 1000). We also show the Hashin-Shtrikman (HS) bounds |3^, the Hori- 
Yonezawa formula which comes from a cumulant series (CS) approximation, and the one-impurity (01) Bruggeman 
formula. The dimension is d = 2. Figs. || and ^ display similar plots for d = 3. 

For moderate contrast (Figs. 1 and 3), we observe that all four self-consistent formulas lie close to each other. This 
is a consequence of the fact that they are exact to second order in the contrast. Also, for any contrast, the slopes 
at p = and p — I are all identical, which is a consequence of the fact that they are exact to second order in the 
dilute limit p — > (the expression near p — 1 is obtained by replacing p hy p — I and by interchanging ai and (T2). 
We also observe that the HS bounds are satisfied in each case considered. However, the formulas obtained via the 
cumulant series summation, i.e. both the HY formula and its "type 2" counterpart, do not reduce to the exact result 
aes = in dimension 1 (not shown). This exact result is also the common value of the HS bounds for d ~ 1. 

Hence, formulas derived from the cumulant series approximation do not obey the HS bounds in dimension d = 1. On 
the other hand, both the Bruggeman formula and its "type 2" counterpart do reduce to the exact result when d = 1, 
and can be seen to always obey the HS bounds whatever d is. The one-impurity approximation scheme therefore 
appears to be of better physical relevance for all dimensions, than the cumulant series approximation. 
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We now discuss the critical behavior. First of all, the percolation thresholds found in the "type 2" formulas are 
Pc = I — exp(— l/(2c?)) (cumulant series) and Pc = l/{2d — 1) (one- impurity). These thresholds are the percolation 
thresholds of the Potts model, and that of the Bethe lattice model, respectively. Both thresholds decrease as l/{2d) 
when d — > oo, which is the exact asymptotics. The Bethe and Potts models are mean-field models, where emphasis 
is put on fluctuations in the couplings between a given site and its neighbours. On the contrary, in effective-medium 
theories, interactions between impurities are taked into account through the self-consistent background medium. Such 
interactions are more important for low dimensions. Effective-medium theories therefore overestimate interactions 
in high dimensions, whereas mean-field models are expected to underestimate them in low dimensions. Above the 
upper critical dimension where mean-field models are accurate, interactions between impurities become irrelevant. 
According to this discussion, our "type 2" formulas appear as hybrids between mean-field and usual effective-medium 
theories, and are expected to be mostly relevant in dimensions intermediate between d — I and the upper critical 
dimension d = 6. Indeed, the condition ai^{ao) = minimizes the influence of the background medium and, according 
to the interpretation developped in Sec. pll B| , replica coupling has to do with couplings between neighboring points. 
The reason for which the introduction of a replica-coupling ansatz yields the exact thresholds of mean-field theories 
will have to be clarified in the future. In Fig. ||, we plot the quadratic fluctuations (Ai?^) / (E)'^ as a function of p. 
The fluctuations in the "type 2" estimates are greatly reduced compared to those of the Bruggeman and HY formula. 
This is consistent with the fact that the influence of the background is reduced. 

"Type 1" formulas give exponents s = t = 1, while for "type 2" formulas they are s = 0, t = 2. Mean-field theories 
yield s = 0, t = 3 which are the exact values for d > 6. It is interesting to compare these values to exact bounds 
deduced from the Nodes-Links-Blobs (NLB) model, in all dimensions. The NLB model is currently accepted as a 
good one for the backbone structure of real random resistor networks The bounds read 

t>l + {d- 2)v, (70a) 
s > 1 + (2 - d)v, (70b) 

where v > Q \s, the correlation length exponent: ^ oc |p — pc\~'^ ■ They hold for 2 < d < 6, whereas for d > 6 the 
right-hand sides in ( [70|) are fixed to their d = 6 values. These bounds follow, e.g., from comparing the lower and upper 
exact bounds obtained in Ref. p3] for the noise exponent n in weakly nonlinear networks, within the NLB scheme. 



They are satisfied by simulation results [g2[. Using the usual effective-medium values s = t = 1 in (70) implies the 
absurd value v — save for d = 2 where a finite value of v is allowed. Though information about the correlation 
length ^ (and therefore about v) is not included in the Bruggeman nor in the HY formulas, the above bounds show 
that, as long as they are meant to model percolating systems obeying the NLB picture, these formulas are truly 
adequate only in dimension d = 2 - and d = 1 where the Bruggeman formula is exact. As to "type 2" formulas, we 
insert the values s = 0, t = 2 into ( [70| ) and deduce that = l/(d — 2), a reasonable expression for d > 3 only. If we 
furthermore insist on having > 1/2 as in real systems, these heuristic arguments restrict the range of validity of the 
new formulas to d = 3,4. Note, moreover, that only in dimension d = 2 are the exponents equal: s = t, because of 
self-duality A formula with unequal exponents therefore is expected to be essentially relevant to dimensions > 3. 

We also quote theoretical bounds for t due to Golden, valid for hierarchical NLB models: 1 < t < 2 for d = 2, 3 
and 2 < t < 3 for d > 4 |3^. The above analysis is consistent with these bounds, and can be summarized as a set of 
prescriptions for using the "best" available effective medium theories, as far as a non-conflicting critical behaviour is 
concerned: for d = 1, Bruggeman's formula, or its "type 2" counterpart are exact; for d = 2, the Bruggeman or HY 
formulas are adequate; for d = 3, 4, "type 2" formulas are applicable; finally, for d > 5 mean-field theories would be the 
most relevant. Estimates or exact values for the exponents are Q: {v,s,t) = (4/3,1.3,1.3)^=2, (0.88,0.73,2.00)^=3, 
(0.68, 0.4, 2.4)d=4, (0.57, 0.1, 2.7)d=5, (1/2, 0, 3)d>6. These values support our prescriptions. 

An interesting observation is that actually both "type 1" and "type 2" formulas can be given by a variational 
formulation as 

a^JJ"' ^ = min crcff(cro), (71) 

<To>0 
0<<!(»0)<1 

cr^yp^ 2 = jnax acff(cro), (72) 

cro>0 
t)<<!(»0)<l 



provided that an unphysical solution (Teff = is discarded in the minimization (71). Indeed, at least in the framework 



of the two different models introduced in Sec. IV, the curves for q{ao) and (Teff(cro) are found to have the form shown 
in Fig. 6. The infimum ( [7l| ) occurs at q = 0, whereas the solution ctq to the equation fT^g(cro) = corresponds to a 
maximum of CTcff. Both types of theories can therefore be interpreted as extremal theories in the framework of self- 
consistent models built on the replica-coupling ansatz. The physical meaning of this interpretation is still not clear. 
However, "type 2" formulas should not been disregarded as unphysical because of their showing up as maximal ones: 
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the minimization principle states that the dissipated power is minimized with respect to the electric field; but there 
is no reason why an extremization with respect to arbitrary variational parameters should not lead to a maximum of 
the dissipated power. Eqs. ([7l|), (72) explain why for a given approximation, one always has cr*g"' ^ < cr*^" ^ in Figs. 
1-4. 

We now consider some points that were not explicitly treated in the paper. First, we presented the formalism in 
terms of the electric field E, from which we obtained a conductivity acs- The electric current j (or the induction D), 
could be used instead [Q. In such a formulation, the random constitutive parameter is the resistivity p{x) — l/a(x) 
and the constraints are V • j = and j — jo. One then computes an effective resistivity pcs- Both formulations are 
equivalent, but a given approximation scheme in general leads to different results for (Teff and aoS =1/ PcB- Preliminary 
investigations of "type 2" formulas have been led in this case. These will be presented elsewhere. Finally, we discuss 
the natural question about the possibility of replica-symmetry breaking [ p5| . Replica symmetry breaking introduces 
more free parameters in the ansatz, and the final extremization has to be carried out with respect to several variables. 
There is no frustration in this problem, and we therefore expect the replica symmetric solution to be the only one. 
In order to test this, we tried a one-step symmetry- breaking solution and did indeed not find any new solution. 



VI. CONCLUSION 



We presented a functional approach to the calculation of effective-medium properties of random media. We showed 
how to recover the Bruggeman and Hori-Yonezawa formulas by using specific approximation schemes to the basic 
functional integral. We also discussed the introduction of a replica-coupling parameter in a gaussian Ansatz, from 
which new effective-medium formulas were obtained. These formulas appear to be more adequate in d = 3 compared 
to the standard ones by Bruggeman and HY. Because it yields a sensible result in all dimensions, and fulfills all 
the constraints required to deserve the label of a "good" effective-medium theory, the "type 2" counterpart of the 
Bruggeman formula offers an interesting alternative to the latter. Indeed, it has a percolation threshold equal to 
Pc = 1/5 in three dimensions. This is closer to values observed in real materials, compared to the pc = 1/3 of the 
Bruggeman formula which often constitutes an overestimation. 
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APPENDIX A: QUADRATIC FLUCTUATIONS OF THE FIELD 

For completeness, we give here the demonstration of Eq. ( |35| ) [ pof . Volume averages are identified to statistical 
ones. Because W* = {l/2){aE'^) — (1/2) YluPv^viE^"^)!^! where denotes an average on the phase v on which the 
conductivity ffi, is constant, the effective conductivity reads 

gcff = ^Py<yy ^^-pl" ■ (Ai) 

On the other hand, deff has to be an homogeneous function of degree one of the CTi,, whence 

cTu-^ — ■ (A2) 

Comparing both equations yields the values of the and consequently that of (i?^) = ^^Pu{E'^)u- Equ. ( |35| ) 

follows. We note that if a^s = (c) (an exact upper bound for the effective conductivity), then (Ai?^) — 0. Therefore, 
Ccff = (c) defines a trivial model of a medium which is a composite, but from which field fluctations are nonetheless 
absent. 
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APPENDIX B: CALCULATIONS IN THE "ONE-IMPURITY" APPROXIMATION 

The approximation which leads to (^) is built as follows. We first expand the exponential 

(^-'''"^-""')„-E^((^-Hof)o. (Bl) 

fe>0 

Since Tio is non-random, using the hamiltonian density wq defined in ( ^8|) and Awx[E{x)] defined by ( ^ ) we can 
rewrite the difference Tie — Ho as 

ne-Ho = J dx AH{x), (B2) 

where 

Ai?(x) = --^ln/e-'^''^"'==[^(^)]\. (B3) 
pv \ I 

The one-impurity approximation consists in writing (fc > 1) 

(He - Uof = J dxi... dxk AH{xi) . . . AH{xk) (B4a) 
~ v'''^ [ dyAH{yf. (B4b) 



The last expression only retains contributions from identical points in Eq. (B4a). Summing back the series in ( |Bl| ), 



and using (B2) yields 

(e-P^n.-n^^))^^l + | ^ _ 1) , (B5) 

f ViE") / '^^ {'^olE{x)]+v Aw^[E{x)]S{x-y)} 

liy) = ^- ^ '- = , (B6) 

where a one-impurity-type integral is involved. Since the fundamental size of the theory v^l'^) is much smaller 
than the volume ^ = 1 of the system, and since the latter is statistically translation- invariant, the outer integral over 
y is redundant with the disorder average, and can be dropped. We therefore arrive at (|o|). 

When /3 oo, the functional can be computed exactly. Let us briefly indicate how to do it. We first 

introduce the notation h for vectors of dimension nd, and components /if, with a — 1, . . . , n, i — 1, . . . , d. Hence, 
Awy\E(^y)\ = Awy{E(yj). The next step is to use the formal identity 

g-/3i,A«;„(S) ^ f -^h■h' -pvAwy{-z-^) ^h'■E /g^N 

J (27r)"'' ^ ' 



to write the numerator of (B6), which we denote hereafter by J{y), as 



(27r)"'' ' 



where M is the matrix defined from the replica-coupling matrix M in wq by M°p = M°'^6ij. After an integration 

over the fields E and (j) implied in the measure 'DE (which can be easily done using the Fourier components of </>, and 
with y — since ( pO[ ) is independent of y), J reads, up to inessential factors [ p9[ : 

Tf,/) = e"^'^/^^^°-^^^°^-^° / '^'^'^'^\ --h-h'^-vPAwy(-,^) ,h'-Eo-h'-AI-^-h'/(20vd) /gg^, 

J (27r)"'^ ■ ^ ' 

Formally expanding exp(—v(3Awy) in powers of —id/dh', and carrying out successive integrations by parts over h' 
yields 
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J{y) 



Det(M) 



'2^ 



d/2 



/ ^fl^~{vl3d/2){Eo-h)-M-{Eo-h)~vf3AiOy{h) 



(BIO) 



where the determinant is evaluated in repHca space. Finally, T{y) = J(y)l J(y \ /S.Wy = 0): 

d/2 



Ay) 



Det(M) 



For any Awy, this integral over the replicated vector field h can be computed exactly using a saddle-point method 
p6t in the limit /3 — > cx), as announced. This allows for a possible extension of the theory to nonlinear media in the 
"one-impurity" approximation. Here, for the linear problem at hand, Eq. (Bll) is a simple gaussian integral. Setting 
Act ~ a — ao and 



Act 



we obtain 



InJ(y) 



-^AafiE^ + ^(In/i - (3vqa^iE^/d) 



(B12) 



(B13) 



from which follows (n2i 



APPENDIX C: CALCULATION IN THE "CUMULANT SERIES" APPROXIMATION 



We have to compute (?ie)o £^nd (7io)o ™ Equ. (jsj). The calculation of (7io)o easy with the methods already 
employed, and yields 



n — >[) Z 



(Cl) 



As to (7ie}o, expand He (Eq. (|T^)) in the cumulants Ck of the disorder averages of cr{x), according to their 

definition by the generating funtion {X is a generic expansion variable) 



fc>i 



(C2) 



We therefore have: 



X k>l 



Cfe(a). 



(C3) 



We deduce that 



^^^)o = -^EE ^r^^^(-)c.(i^V2), 

^ X A;>1 



(C4) 



where 



Cfc(i?V2) 



E^"(^)'/2 



(C5) 



(because of statistical homogeneity, these coefficients do not depend on the position variable x). It is convenient to 
introduce the following generating function Z{X) in order to compute the Ck- 
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k>l 



1 

n 



cxp 



- 1 



— In ( exp 
n \ 



0[n) 



(C6) 



Setting = + (X/z;/3d)[M-i]"T, we obtain [p9| 



Z(X) = J dTihnA + XElY^[A-^]'^'' \ 

K a-y ) 



0{n) 



(C7) 



(the trace and the logarithm act in the rephca space). 

Expanding (C7) in powers of X then aUows for the identification 





{— 








1 n 



a-y 



Oin). 



Use of this expression in (C4) cancels the convergence factor kl: we reintroduce it by inserting the identity 



1 



I due-^'u"' = 1, 

r' 



(C8) 



(C9) 



applied to m fc — 1 and to = fc in the resulting cumulant series. This permits its Borel summation, which brings in 
the functions hm{x) defined by (Cll). This results in 



1 



1 



ElY,[Mh,{M-^/d)r + ^Tr /i_i(Af-Vd) 



where we defined the family of functions 



hm{z) 



duu'"e-"ln(e-'""^) (to, > -2). 



0(n) 



Note for further use that 



h'mi^) = - [hm+liz) - (fc + l)hrniz) 
Z 



(CIO) 



(Cll) 



(C12) 



For M given by (^9|), the differents terms in (|C10| ) are 

lim - Y^[Mho{M-^/d)Y^ = ao/io(l/(rf'To)), 



a7 



lim -TyK.AM-^ Id) = dh^AXlidao)) + v(i da^qhMl {daQ))El. 



(C13) 



which leads to Eq. (p3 
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FIG. 1. Rescaled effective conductivities in dimension 
d = 2 for a binary medium, versus tfie volume concentation 
p of component 2. Tlie conductivity ratio is a^jai = 10. 
Higfiest and lowest solid curves: Hashin-Shtrikman bounds; 
Br.: the Bruggeman formula ("type 1", one-impurity approx- 
imation - 01); HY: the Hori-Yonezawa formula ("type 1", 
cumulant series approximation - CS); both "type 2" curves 
are the new formulas, within OI and CS approximations. 
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FIG. 3. Rescaled effective conductivities in dimension 
d = 3 for a binary medium, versus the volume concentation p 
of component 2. The conductivity ratio is a^jox — 10. Same 
plots as in Fig. 0. 




FIG. 2. Rescaled effective conductivities in dimension 
d = 2 for a binary medium, versus the volume concentation 
p of component 2. The conductivity ratio is a-ijax = 1000. 
Same plots as in Fig. hi 



FIG. 4. Rescaled effective conductivities in dimension 
d = 3 for a binary medium, versus the volume concentation 
p of component 2. The conductivity ratio is a^jax ~ 1000. 
Same plots as in Fig. m 
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FIG. 5. Relative quadratic fluctuations of the fleld in di- 
mension d — 3 for a binary medium, versus the volume con- 
centation p of component 2. The conductivities are ai — 1, 
(72 = 1000. 
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FIG. 6. Rescaled effective conductivity (Tcff((7o)/(Ti, and re- 
duced replica coupling parameter g(cro) vs. uq in dimension 
d — 3 for a binary medium. In this exemple computed from 
Eqs. (^^, the conductivity ratio is a2/ui = 100, and 

the volume fraction p of component 2 is p = 0.18. The "type 
1" effective conductivity is obtained when q = 0, whereas the 
"type 2" effective conductivity corresponds to the maximum 
of the curve aefi(o"o). 
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